Transcriptome Analysis of Tomato Leaves Reveals Candidate Genes Responsive to Tomato Brown Rugose Fruit Virus Infection.

Tomato brown rugose fruit virus (ToBRFV) is a newly-emerging tobamovirus which was first reported on tomatoes in Israel and Jordan, and which has now spread rapidly in Asia, Europe, North America, and Africa. ToBRFV can overcome the resistance to other tobamoviruses conferred by tomato Tm-1, Tm-2, and Tm-22 genes, and it has seriously affected global crop production. The rapid and comprehensive transcription reprogramming of host plant cells is the key to resisting virus attack, but there have been no studies of the transcriptome changes induced by ToBRFV in tomatoes. Here, we made a comparative transcriptome analysis between tomato leaves infected with ToBRFV for 21 days and those mock-inoculated as controls. A total of 522 differentially expressed genes were identified after ToBRFV infection, of which 270 were up-regulated and 252 were down-regulated. Functional analysis showed that DEGs were involved in biological processes such as response to wounding, response to stress, protein folding, and defense response. Ten DEGs were selected and verified by qRT-PCR, confirming the reliability of the high-throughput sequencing data. These results provide candidate genes or signal pathways for the response of tomato leaves to ToBRFV infection.

First report of citrus leaf blotch virus infecting Forsythia viridissima in China.

Green-stem forsythia (Forsythia viridissima), also known as golden bell, is cultivated widely in China as an early spring flowering shrub. In July 2020, yellow or white vein clearing symptoms on leaves were observed in approximate 15% golden bell plants along a landscape river in Ningbo city, Zhejiang province, China. Symptomatic leaves from six different plants were collected and pooled. Total RNA was extracted from about 200 mg pooled sample using TRIzol Reagent (Invitrogen, Carlsbad, USA) and used for high-throughput sequencing (HTS). The cDNA library was constructed using a TruSeq RNA Sample Preparation Kit (Illumina) and an Illumina NovaSeq 6000 platform was utilized to yield 150 nt paired-end reads. CLC Genomic Workbench 11 (QIAGEN) with default parameters were used for data analysis. A total of 41,604,174 paired-end reads were obtained, and 156,853 contigs (16 - 26,665 nt) were generated de novo and compared with sequences in the NCBI nt and nr database using BLASTn and BLASTx, respectively. A total of 197,277 reads were mapped to the citrus leaf blotch virus (CLBV; genus Citrivirus, family Betaflexiviridae) genome with an average coverage of 3191×. A contig of 8783 nt (excluding the poly(A) tail) was aligned to CLBV isolate Vib (accession No. OP751940) by BLASTn with the highest nt sequence identity of 99.7% and 99% query coverage, suggesting that the samples were infected with CLBV (Myung-Hwi Kim et al. 2023). No other virus was detected by this analysis. Subsequently, leaves of the six plants collected above, three plants with mild chlorotic symptoms and three plants without obvious symptoms were tested separately by RT-PCR and all were positive for CLBV. Sap from multiple symptomatic F. viridissima leaves was mechanically inoculated to Nicotiana benthamiana, N. tabacum and Datura stramonium in sextuplicate, but after two months, none of the inoculated plants had obvious symptoms and all of them tested negative for CLBV using RT-PCR. To determine the genome sequence of CLBV present in F. viridissima, a single sample from one plant was selected for genome validtion. The contig sequence was confirmed by Sanger sequencing of RT-PCR products amplified using CLBV-specific primers, and the 5' terminal sequence of the virus was determined using a commercial SUPERSWITCH RACE cDNA Synthesis Kit (Tiosbio, Beijing, China). The complete genomic sequence of CLBV isolated from F. viridissima was 8787 nts long, excluding the poly(A) tail, has the expected three predicted ORFs and was deposited in the GenBank database (accession no. OR766026). Phylogenetic analysis of different CLBV genome sequences from fruit trees and other hosts in GenBank using MEGA11 showed that the golden bell isolate was most closely related to isolate Vib (OP751940) from Viburnum lentago in South Korea, with which it was almost identical (99.7% complete nt sequence identity and >99% aa sequence identity in each of the three ORFs). Ten viruses have been previously reported from Forsythia spp. (Kaminska, M. 1985; Lee et al. 1997), but this is the first report of CLBV in this host. CLBV mainly infects citrus, kiwifruit and apple causing mosaic, chlorosis or yellow vein clearing symptoms, however, bud union disorder was observed in 'Nagami' kumquat infected by CLBV, which caused serious production losses (Cao et al. 2017; Li et al. 2018; Liu et al. 2019; Galipienso et al. 2001). Therefore, further investigation is needed to assess if F. viridissima can be an intermediate host to transfer CLBV to other crops.

Adoption of the 2A Ribosomal Skip Principle to Track Assembled Virions of Pepper Mild Mottle Virus in Nicotiana benthamiana.

The coat protein (CP) is an important structural protein that plays many functional roles during the viral cycle. In this study, the CP of pepper mild mottle virus (PMMoV) was genetically fused to GFP using the foot-and-mouth disease virus peptide 2A linker peptide and the construct (PMMoV-GFP2A) was shown to be infectious. The systemic spread of the virus was monitored by its fluorescence in infected plants. Electron microscopy and immunocolloidal gold labelling confirmed that PMMoV-GFP2A forms rod-shaped particles on which GFP is displayed. Studies of tissue ultrastructure and virion self-assembly confirmed that PMMoV-GFP2A could be used to monitor the real-time dynamic changes of CP location during virus infection. Aggregations of GFP-tagged virions appeared as fluorescent plaques in confocal laser microscopy. Altogether, PMMoV-GFP2A is a useful tool for studying the spatial and temporal changes of PMMoV CP during viral infection.

Integrated ATAC-seq and RNA-seq data analysis identifies transcription factors related to rice stripe virus infection in Oryza sativa.

Animal studies have shown that virus infection causes changes in host chromatin accessibility, but little is known about changes in chromatin accessibility of plants infected by viruses and its potential impact. Here, rice infected by rice stripe virus (RSV) was used to investigate virus-induced changes in chromatin accessibility. Our analysis identified a total of 6462 open- and 3587 closed-differentially accessible chromatin regions (DACRs) in rice under RSV infection by ATAC-seq. Additionally, by integrating ATAC-seq and RNA-seq, 349 up-regulated genes in open-DACRs and 126 down-regulated genes in closed-DACRs were identified, of which 34 transcription factors (TFs) were further identified by search of upstream motifs. Transcription levels of eight of these TFs were validated by reverse transcription-PCR. Importantly, four of these TFs (OsWRKY77, OsWRKY28, OsZFP12 and OsERF91) interacted with RSV proteins and are therefore predicted to play important roles in RSV infection. This is the first application of ATAC-seq and RNA-seq techniques to analyse changes in rice chromatin accessibility caused by RSV infection. Integrating ATAC-seq and RNA-seq provides a new approach to select candidate TFs in response to virus infection.

The nanovirus U2 protein suppresses RNA silencing via three conserved cysteine residues.

Nanoviruses have multipartite, circular, single-stranded DNA genomes and cause huge production losses in legumes and other crops. No viral suppressor of RNA silencing (VSR) has yet been reported from a member of the genus Nanovirus. Here, we demonstrate that the nanovirus U2 protein is a VSR. The U2 protein of milk vetch dwarf virus (MDV) suppressed the silencing of the green fluorescent protein (GFP) gene induced by single-stranded and double-stranded RNA, and the systemic spread of the GFP silencing signal. An electrophoretic mobility shift assay showed that the U2 protein was able to bind double-stranded 21-nucleotide small interfering RNA (siRNA). The cysteine residues at positions 43, 79 and 82 in the MDV U2 protein are critical to its nuclear localization, self-interaction and siRNA-binding ability, and were essential for its VSR activity. In addition, expression of the U2 protein via a potato virus X vector induced more severe necrosis symptoms in Nicotiana benthamiana leaves. The U2 proteins of other nanoviruses also acted as VSRs, and the three conserved cysteine residues were indispensable for their VSR activity.

Complete genome sequence of Paris polyphylla chlorotic mottle virus infecting Paris polyphylla var. yunnanensis in southwest China.

A novel virus infecting a Paris polyphylla var. yunnanensis plant, tentatively named "Paris polyphylla chlorotic mottle virus" (PpCMV), was discovered in the city of Lijiang, Yunnan Province, China. Its genome consists of 6384 nucleotides (nt), excluding the 3'-terminal poly(A) tail, and contains two open reading frames: ORF1 and ORF2. ORF1 is 6150 nt in length, encoding a large 2050-aa polyprotein with at least two conserved regions encoding a replication-associated protein and a coat protein, the latter of which is located at the 3' end of ORF1. ORF2, consisting of 1185 nt, is located within ORF1 but has a different reading frame. It encodes a 394-aa-long putative movement protein. Phylogenetic analysis based on amino acid sequences revealed that the newly discovered virus exhibited the closest relationship to Hobart betaflexivirus 1 and rhodiola betaflexivirus 1, both of which belong to the genus Capillovirus, sharing 48.8% and 36.5% amino acid sequence identity, respectively, in the structural protein. This is the first report of the complete genome sequence of PpCMV in China.

Complete genome sequence of a novel mitovirus detected in Colocasia esculenta.

A novel mitovirus was detected in taro (Colocasia esculenta) growing in Ningbo, China. The complete genome sequence of Colocasia esculenta associated mitovirus 1 (CeaMV1) was determined by next-generation sequencing combined with RT-PCR and RACE. The genome is 2921 nucleotides long and contains a single ORF encoding a putative RNA-dependent RNA polymerase. Homology searches and phylogenetic analysis suggested that CeaMV1 is a member of a new species in the genus Duamitovirus. This is the first report of a member of the family Mitoviridae associated with taro.

Pepper mild mottle virus can infect and traffick within Nicotiana benthamiana plants in non-virion forms.

Virions are responsible for the long-distance transport of many viruses, such as Pepper mild mottle virus (PMMoV). Emerging evidence indicates viral traffic in the form of ribonucleoprotein complexes (RNP), yet comprehensive analysis is scarce. In this study, we inoculated plants with PMMoV-GFP, both with and without the coding sequence for the coat protein (CP). PMMoV-GFP was detected in systemic leaves, even in the absence of the CP, despite the presence of much smaller infection areas. Moreover, using leaf extracts from PMMoV-infected plants to perform a root-irrigation experiment, we confirmed that PMMoV can infect plants through root transmission. Diluting the leaf extracts significantly diminished infectivity, and attempts to compensate for the dilution of other components by adding virions above the original level proved ineffective. Our findings strongly indicate that PMMoV can infect and traffick within plants in non-virion forms. Future studies should aim to identify the specific forms involved.

The Additional 15 nt of 5' UTR in a Novel Recombinant Isolate of Chilli Veinal Mottle Virus in Solanum nigrum L. Is Crucial for Infection.

An isolate of chilli veinal mottle virus (ChiVMV; genus Potyvirus) of Solanum nigrum L. from southwest China (ChiVMV-YunN/Yuxi) was identified and sequenced (GenBank: OP404087). Comparison with other ChiVMV isolates and recombination analyses suggested a recombinant origin. The most significant recombination event among all 21 complete ChiVMV isolates was an ending breakpoint at 1408-1488 for ChiVMV-YunN/Yuxi with ChiVMV-TaiW and ChiVMV-YunN/Ca operating as the respective major and minor parents. Interestingly, the 5' UTR of ChiVMV-YunN/Yuxi is 15 nucleotides ('AAAAATAAAACAACC') longer than other reported isolates. A full-length clone of ChiVMV-YunN/Yuxi was constructed and was shown to be infectious in Nicotiana benthamiana. The additional 15 nt of 5' UTR in ChiVMV-YunN/Yuxi was stable when transmitted through three generations. Experiments with modified clones showed that the additional 15 nt are essential for infection by this isolate.

Selection and Validation of Reference Genes in Virus-Infected Sweet Potato Plants.

Quantitative real-time PCR (qRT-PCR) in sweet potatoes requires accurate data normalization; however, there are insufficient studies on appropriate reference genes for gene expression analysis. We examined variations in the expression of eight candidate reference genes in the leaf and root tissues of sweet potatoes (eight nonvirus-infected or eight virus-infected samples). Parallel analyses with geNorm, NormFinder, and Best-Keeper show that different viral infections and origin tissues affect the expression levels of these genes. Based on the results of the evaluation of the three software, the adenosine diphosphate-ribosylation factor is suitable for nonvirus or virus-infected sweet potato leaves. Cyclophilin and ubiquitin extension proteins are suitable for nonvirus-infected sweet potato leaves. Phospholipase D1 alpha is suitable for virus-infected sweet potato leaves. Actin is suitable for roots of nonvirus-infected sweet potatoes. Glyceraldehyde-3-phosphate dehydrogenase is suitable for virus-infected sweet potato roots. The research provides appropriate reference genes for further analysis in leaf and root samples of viruses in sweet potatoes.

Complete genome sequence of a new virus from Allium sativum L in China.

The complete genome of a new virus belonging to the family Betaflexiviridae was identified in garlic and sequenced by next-generation sequencing and reverse transcription PCR. The complete RNA genome (GenBank accession number OP021693) is 8191 nucleotides in length, excluding the 3' poly(A) tail, and contains five open reading frames (ORFs). These open reading frames encode the viral replicase, triple gene block, and coat protein, and the genome organization is typical of members of the subfamily Quinvirinae. The virus has been tentatively named "garlic yellow curl virus" (GYCV). Phylogenetic analysis suggested that it represents an independent evolutionary lineage in the subfamily, clustering with the currently unclassified garlic yellow mosaic associated virus (GYMaV) and peony betaflexivirus 1 (PeV1). Differences between the phylogenies inferred for the replicase and coat protein indicate that the new virus does not belong to any established genus of the family Betaflexiviridae. This is the first report of GYCV in China.

Pepper mild mottle virus coat protein interacts with pepper chloroplast outer envelope membrane protein OMP24 to inhibit antiviral immunity in plants.

Pepper mild mottle virus (PMMoV) is a devastating viral pathogen of pepper (Capsicum annuum) but it is unclear whether and how peppers protect against PMMoV infection. The expression of the chloroplast outer membrane protein 24 (OMP24) of C. annuum was upregulated under PMMoV infection and it interacted with PMMoV coat protein (CP). Silencing of OMP24 in either C. annuum or Nicotiana benthamiana facilitated PMMoV infection, whereas overexpression of N. benthamiana OMP24 in transgenic plants inhibited PMMoV infection. Both C. annuum OMP24 (CaOMP24) and N. benthamiana OMP24 (NbOMP24) localized to the chloroplast and have a moderately hydrophobic transmembrane domain that is necessary for their localization. Overexpression of CaOMP24 induced stromules, perinuclear chloroplast clustering, and accumulation of reactive oxygen species (ROS), the typical defense responses of chloroplasts transferring the retrograde signaling to the nucleus to regulate resistance genes. The expression of PR1 and PR2 was also upregulated significantly in plants overexpressing OMP24. Self-interaction of OMP24 was demonstrated and was required for OMP24-mediated plant defense. Interaction with PMMoV CP interfered with the self-interaction of OMP24 and impaired OMP24-induced stromules, perinuclear chloroplast clustering and ROS accumulation. The results demonstrate the defense function of OMP24 in pepper during viral infection and suggest a possible mechanism by which PMMoV CP modulates the plant defense to facilitate viral infection.

Complete nucleotide sequence of hackberry virus A, a tentative member of the genus Waikavirus.

The complete genomic sequence of a waikavirus from Chinese hackberry in Zhejiang province, China, named "hackberry virus A" (HVA), was determined using high-throughput sequencing (HTS) combined with reverse transcription polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE) PCR. The bicistronic genomic RNA of HVA was found to consist of 12,691 nucleotides (nt), excluding the 3'-terminal poly(A) tail, and to encode a large polyprotein of 3783 amino acids (aa) and an additional 10.3-kDa protein. The aa sequences of the Pro-Pol and the CP regions of this virus share 39.8-44.2% and 25.5-36.4% identity, respectively, with currently known waikaviruses. These values are significantly below the current species demarcation threshold (< 75% and < 80% aa identity for the CP and Pro-Pol region, respectively) for the family Secoviridae, indicating that HVA represents a new species in the genus Waikavirus. This is the first report of a virus infecting Chinese hackberry.

Selection and Validation of Reference Genes for RT-qPCR Analysis of Gene Expression in Nicotiana benthamiana upon Single Infections by 11 Positive-Sense Single-Stranded RNA Viruses from Four Genera.

Quantitative real-time PCR (RT-qPCR) is a widely used method for studying alterations in gene expression upon infections caused by diverse pathogens such as viruses. Positive-sense single-stranded (ss(+)) RNA viruses form a major part of all known plant viruses, and some of them are damaging pathogens of agriculturally important crops. Analysis of gene expression following infection by ss(+) RNA viruses is crucial for the identification of potential anti-viral factors. However, viral infections are known to globally affect gene expression and therefore selection and validation of reference genes for RT-qPCR is particularly important. In this study, the expression of commonly used reference genes for RT-qPCR was studied in Nicotiana benthamiana following single infection by 11 ss(+) RNA viruses, including five tobamoviruses, four potyviruses, one potexvirus and one polerovirus. Stability of gene expression was analyzed in parallel by four commonly used algorithms: geNorm, NormFinder, BestKeeper, and Delta CT, and RefFinder was finally used to summarize all the data. The most stably expressed reference genes differed significantly among the viruses, even when those viruses were from the same genus. Our study highlights the importance of the selection and validation of reference genes upon different viral infections.

Complete genome sequence of iris potyvirus B infecting Lilium lancifolium in China.

The complete genome sequence of a virus from lily (Lilium lancifolium Thunb.) growing in Huoshan County, Anhui Province, China, was determined. The whole genome consists of 9558 nucleotides, excluding the poly(A) tail, and encodes a 3061-amino-acid polyprotein (GenBank number ON365558) typical of potyviruses. This is the first complete genome sequence of iris potyvirus B (IPB), for which only a partial sequence from Iris domestica was reported previously. Comparative analysis of this genome sequence with those of closely related potyviruses identified nine cleavage sites and the conserved motifs typical of potyviruses. The complete polyprotein ORF shares 73.6% nucleotide and 81.6% amino acid sequence identity with that of iris potyvirus A (IPA, GenBank number MH898493). Phylogenetic analysis showed that IPB is related to IPA and clusters in a group with lily yellow mosaic virus (LYMV). This is the first report of IPB infecting lily plants.

Construction of an infectious full-length and eGFP-tagged cDNA clone of a chilli ringspot virus isolate from Yunnan province, China.

Chilli ringspot virus (ChiRSV; genus Potyvirus) was one of several viruses previously detected in pepper samples with severe yellowing and curling symptoms growing in Wenshan, Yunan province, China. We now report the full-length sequence of ChiRSV-YN/Wenshan (MZ269480), which has 88.5-98.9% nucleotide sequence identity to other published ChiRSV isolates. A full-length cDNA infectious clone was constructed. This cDNA and an eGFP-tagged clone were infectious, leading to systemic symptoms in both Nicotiana benthamiana and Capsicum spp. Recombinant clones containing the P1 protein coding region of other ChiRSV isolates differed in their pathogenicity. Single infection by ChiRSV caused mild mosaic or leaf crinkling in Capsicum frutescens L. and Capsicum annuum L.

Acidic dileucine motifs in the cylindrical inclusion protein of turnip mosaic virus are crucial for endosomal targeting and viral replication.

Previously we reported that the multifunctional cylindrical inclusion (CI) protein of turnip mosaic virus (TuMV) is targeted to endosomes through the interaction with the medium subunit of adaptor protein complex 2 (AP2ß), which is essential for viral infection. Although several functionally important regions in the CI have been identified, little is known about the determinant(s) for endosomal trafficking. The CI protein contains seven conserved acidic dileucine motifs [(D/E)XXXL(L/I)] typical of endocytic sorting signals recognized by AP2ß. Here, we selected five motifs for further study and identified that they all were located in the regions of CI interacting with AP2ß. Coimmunoprecipitation assays revealed that alanine substitutions in the each of these acidic dileucine motifs decreased binding with AP2ß. Moreover, these CI mutants also showed decreased accumulation of punctate bodies, which enter endocytic-tracking styryl-stained endosomes. The mutations were then introduced into a full-length infectious clone of TuMV, and each mutant had reduced viral replication and systemic infection. The data suggest that the acidic dileucine motifs in CI are indispensable for interacting with AP2ß for efficient viral replication. This study provides new insights into the role of endocytic sorting motifs in the intracellular movement of viral proteins for replication.

Turnip mosaic virus co-opts the vacuolar sorting receptor VSR4 to promote viral genome replication in plants by targeting viral replication vesicles to the endosome.

Accumulated experimental evidence has shown that viruses recruit the host intracellular machinery to establish infection. It has recently been shown that the potyvirus Turnip mosaic virus (TuMV) transits through the late endosome (LE) for viral genome replication, but it is still largely unknown how the viral replication vesicles labelled by the TuMV membrane protein 6K2 target LE. To further understand the underlying mechanism, we studied the involvement of the vacuolar sorting receptor (VSR) family proteins from Arabidopsis in this process. We now report the identification of VSR4 as a new host factor required for TuMV infection. VSR4 interacted specifically with TuMV 6K2 and was required for targeting of 6K2 to enlarged LE. Following overexpression of VSR4 or its recycling-defective mutant that accumulates in the early endosome (EE), 6K2 did not employ the conventional VSR-mediated EE to LE pathway, but targeted enlarged LE directly from cis-Golgi and viral replication was enhanced. In addition, VSR4 can be N-glycosylated and this is required for its stability and for monitoring 6K2 trafficking to enlarged LE. A non-glycosylated VSR4 mutant enhanced the dissociation of 6K2 from cis-Golgi, leading to the formation of punctate bodies that targeted enlarged LE and to more robust viral replication than with glycosylated VSR4. Finally, TuMV hijacks N-glycosylated VSR4 and protects VSR4 from degradation via the autophagy pathway to assist infection. Taken together, our results have identified a host factor VSR4 required for viral replication vesicles to target endosomes for optimal viral infection and shed new light on the role of N-glycosylation of a host factor in regulating viral infection.